A precise method for visualizing dispersive features in image plots 
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In order to improve the advantages and the reliability of the second derivative method in tracking 
the position of extrema from experimental curves, we develop a novel analysis method based on 
the mathematical concept of curvature. We derive the formulas for the curvature in one and two 
dimensions and demonstrate their applicability to simulated and experimental angle-resolved pho- 
toemission spectroscopy data. As compared to the second derivative, our new method improves the 
localization of the extrema and reduces the peak broadness for a better visualization on intensity 
image plots. 
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I. INTRODUCTION 

With the development of multi-channel detectors and 
the recording of a huge amount of experimental data, 
the pass decade has witnessed a boom in the use of color 
images for the representation of spectroscopic data in a 
very compact and easily visualized way. Typically, a color 
scale is associated with the experimental spectral inten- 
sity, which is displayed as a function of two independent 
variables. For example, such images are widely used 
in scanning^ tunneling microscopy (STM) Raman 
scattering 0-01 > inelastic neutron scattering (INS) @-[ll[ , 
atomic force microscopy (AF M) [12HL4j , resonant inelas- 
tic X-ray scattering (RIXS) jl5l - fl8| and ang le-resolved 
photoemission spectroscopy (ARPES) |19l - [30| . 

This imaging process is particularly efficient to rep- 
resent energy band dispersions in the momentum or 
momentum-transfer spaces, where the energy and the 
momentum (or momentum-transfer) are the two inde- 
pendent variables. Frequently though, many bands or 
features overlap or have significant broadness, making di- 
rect visualization of the raw data difficult. The main tool 
commonly used in ARPES analysis to overcome this issue 
and improve direct visualization of band disp ersion is the 
second derivative of intensity plots |20l - (29j . Despite its 
success and widespread use, the method of second deriva- 
tive gives sometimes results that differs slightly from the 
actual position of the maxima in the energy distribu- 
tion curves (EDCs), where the photoemission intensity at 
fixed momentum is represented as a function of energy, 
or in the momentum distribution curves (MDCs), where 
the photoemission intensity at fixed energy is given as a 
function of momentum. Alternatives must thus be found 
to improve both accuracy and visualization of data. 

In this paper, we develop an analysis method for 
studying spectroscopic data based on the mathemati- 
cal concept of curvature in one-dimension (ID) and two- 
dimension (2D). As an example, we apply this method 
to the study of electronic energy dispersion from ARPES 
data. We show two major advantages of the curvature 
method over the second derivative method: (i) the cur- 



vature method is more reliable in tracking the position of 
extrema and (ii) the curvature method can increase the 
sharpness of the dispersive features for a better visualiza- 
tion effect. We prove the efficiency of this method using 
both experimental and simulated data. 



II. ID CURVATURE METHOD 

The concept of curvature is used to quantitatively de- 
termine how much a curve is not straight. It locally asso- 
ciates a radius of curvature, which can be either positive 
or negative, to a small segment along a curve. The math- 
ematical definition in ID of the curvature C(x) associated 
to a function f(x) is given by: 
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For application to spectroscopic data, for example to 
an EDC curve, f(x) may represent the signal intensity 
whereas x represents a unitlcss variable such as normal- 
ized energy. The normalization of a variable x that 
carries units is done through a transformation such as 
x/£ — > x, where £ is a positive arbitrary constant with 
the same dimension as x. Since experimental spectro- 
scopic functions themselves are usually defined to an arbi- 
trary factor, f(x) carries the same information as I f(x), 
where Iq is an arbitrary positive constant. Taking into 
account the arbitrariness in the absolute values of x and 
f(x), we can rewrite equation ([1]) as: 
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Since we are interested uniquely in the relative variations 
of the curvature, this equation can be reduced further to: 
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FIG. 1. Color online, (a) Simulated ARPES intensity plot (see the text), (b) Spectral intensity (red curve) as a function of 
momentum (MDC) along the horizontal red line in panel (a), compared to the curvature with different values of ao (see the 
text), (c) Comparison between the spectral intensity as a function of energy (EDC) along the vertical black line in panel (a), 
and the corresponding second derivative and curvature curves, (d) Same as (c) but for the spectral intensity as a function of 
momentum along the red line in panel (a). The intensity of each curve in panels (b)-(d) has been normalized to 1 and the sign 
of the second derivative and curvature curves has been reversed to facilitate visualization. The intensity plots of the second 
derivatives [curvature] of the simulated data from panel (a) along the energy and momentum directions are given in (e) [(g)] 
and (f) [(h)], respectively. 



where Co is a free parameter. In order to understand the 
meaning of Co , we test the previous equation in two limit 
cases: 



(1) When Co > f(x) 



/hen f'(x) 2 can be ig- 



nored, we get 
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which gives the same result as the second derivative 
method. 

(2) When Co < f'{xf, 
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This latter solution diverges at the extrema, where 
f'{x) = 0. As Co approaches 0, the peak positions in 
C{x) are getting closer and closer to the real peak posi- 
tions. In the worst case, when Co — > oo, the curvature 
should provide a result as good as the one given by the 
second derivative. Therefore, the curvature is necessar- 
ily an improvement over the second derivative method in 



tracking the peak positions. In practice, we avoid sin- 
gularities while maintaining the reliability of C(x) by 
choosing an intermediate Co- Empirically, we find out 
that the best compromise is reached when Co is of the 
order of the average or the maximum value of |/'(a;)| 2 . 
Hereafter, we express Co as ao|/'(a;) | \ lax - l where ao is a 
positive constant and \f \x)\ ma x is the maximum value 
of \f(x)\. 

To illustrate the reliability of the curvature analysis, 
we simulate ARPES data using known parameters. The 
ARPES photoemission intensity can be expressed by the 
product of three terms: the Fermi-Dirac distribution 
Jd{x), the spectral weight A(k, oj) that contains all the 
information about the dispersion, and a matrix element 
factor that depends on momentum, as well as on the 
energy and polarization of the probing photons. Since 
the latter term does not carry any information about the 
dispersion, we set it to 1. The spectral weight can be 
expressed in terms of the energy dispersion £k as: 
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where S(k, cj)=E'(k, u;)+i£"(k, oj) is the self-energy of 
the quasi-particles. The self-energy is known to depend 
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only weakly on momentum and its imaginary part usually 



varies like < 
energy to: 



-c at low energy. Thus, we set the self- 
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which satisfies the Kramers-Kronig transformation [3l| . 

Setting a = 3 and c = 0.15 eV, we plot simulated 
ARPES data in Fig. [Ha) for the dispersion eu = 
15k 2 - 0.3 eV at a temperature (T) of 20 K. The result 
has been further convoluted by a Gaussian function along 
the energy direction to simulate an energy resolution of 
10 meV. In Fig. QJ b), we compare the MDC along the red 
line in panel (a) to curvature curves of that same MDC 
using different values of ao . For a better comparison, the 
sign of the curvature curves has been reversed and the 
maxima of all curves have been normalized to 1. As ex- 
pected for an asymmetrical lineshape, the position of the 
curvature peak is slightly away from the real peak posi- 
tion when ao is large but converges to that latter position 
with ao decreasing. Moreover, the peak sharpens rapidly 
as a decreases. Although this is obviously an advantage 
in tracking its position, we note that it is necessary to re- 
frain decreasing ao too much while studying multi-feature 
systems since the sharpening of the peaks is accompanied 
by an increase of intensity in the curvature, which may 
affect the global contrast between all the features repre- 
sented on a single image. We also note that since we are 
trying to find peak positions (maxima or inflections in 
the spectra), only the positive parts of the sign-reversed 
second derivatives and the sign-reversed curvatures have 
a physical meaning (the approximate position of peaks), 
and the negative parts are completely ignored. 

In Fig. rile) and Fig.ETJd), we plot the EDC and MDC 
along the black and red lines in panel (a), respectively, 
along with their second derivative and curvature curves 
(normalized and sign- reversed). Since both the MDC 
and EDC lineshapes are asymmetric with respect to the 
peak positions, the second derivative curves do not track 
the peak positions exactly and a small shift towards the 
highest slope change is observed. In contrast, the cur- 
vature analysis provides more reliable peak positions, in 
addition to giving sharper features. We performed the 
second derivative analysis for all EDCs and MDCs and 
we show the corresponding second derivative intensity 
plots in Fig. [TJe) and Fig. [TJf ) , respectively. Similarly, 
the EDC- and MDC-curvature intensity plots associated 
with the data of Fig. [IJa) are given in Fig. [IJg) and 
Fig.QJh), respectively. Obviously, the curvature method 
gives sharper features and allows a better tracking of the 
band dispersion as compared with the second derivative 
analysis. However, as for the analysis of EDCs and MDCs 
and their corresponding second derivatives, the ID cur- 
vature method presented here has some limitations over 



the whole range of energy and momentum. While the 
EDC-curvature method is quite reliable to track the min- 
ima and maxima of band dispersions, it gives unreliable 
results near the Fermi cutoff, which itself appears as a 
spectral feature. In contrast, the MDC-curvature method 
is quite precise near the Fermi cutoff but fails to reveal 
precisely the dispersion near extrema. Nevertheless, a 
cleaver combined use of EDC- and MDC-curvature anal- 
ysis allows to track the band dispersion completely and 
precisely. A more sophisticated analysis method is pro- 
posed in the next section. 

We now test the ID curvature method on real exper- 
imental data. In Fig. [U[a), we show an intensity plot 
recorded at 15 K corresponding to the low-energy band 
dispersion near the Fermi wavevector (kp) of the so- 
called a band in optimally-doped Bao.6Ko.4Fc2As2 (T c 
= 37 K) [lH. As reported earlier, the dispersion ex- 
hibits in the superconducting state a kink or sudden 
slope change around 25 meV below the Fermi energy 
(Ep) due to an electron- mode coupling [33}. Although 
the kink is visible in the original image, it appears more 
clearly in the MDC-second derivative plot shown in Fig. 
[2{b). As expected from the previous discussion, the re- 
sult is even sharper with the use of the MDC-curvature 
method, as illustrated in Fig. [U(c). The second derivative 
method is particularly efficient in ARPES for the study 
of band dispersion complexes. In Fig. HJd), we show an 
ARPES intensity cut of Sr4V20 6 Fe 2 As2 recorded at 40 
K along the T— M direction J34|. Within the wide energy 
range displayed (down to about 1.5 eV below Ep), many 
bands exist and overlap, and it is very difficult to extract 
their band dispersion. The corresponding EDC-second 
derivative intensity plot shown in Fig. [2Ie) is a clear 
improvement for the visualization of the main bands. 
Once more, this advantage is reinforced with the EDC- 
curvature method, as illustrated in Fig.[5ff). The bands 
are sharper and the reliability in tracking the peak posi- 
tion improved. 



III. 2D CURVATURE METHOD 

Despite its ability to track band dispersions, the ID 
curvature method has some unavoidable problems when 
analyzing intensity images. The main problem comes 
from the fact that the images themselves, as well as the 
features they emphasize, are 2D rather than ID objects. 
In this section, we extend the ID curvature method to a 
2D method. As a first example, we treat the simplified 
case where the two independent variables determining 
the spectral intensity are equivalent. For example, this 
situation applies to AFM and STM mappings, for which 
both independent variables represent a distance, as well 
as to ARPES Fermi surface mappings, for which both in- 
dependent variables represent a momentum component. 
Aftcwards, we will focus on a more general case, where 
the independent variables are inequivalent, like in the 
energy vs momentum intensity plots used in ARPES to 
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FIG. 2. Color online, (a) ARPES intensity plot (from [331). 
(b)[(c)] Corresponding intensity plot of second derivative [ID 
curvature] alon g th e momentum direction, (d) ARPES inten- 
sity plot (from [331). (e)[(f)] Corresponding intensity plot of 
second derivative [ID curvature] along the energy direction. 



reveal energy band dispersions. 



A. Equivalent independent variables 

The equivalent in 2D of the second derivative is the 
Laplacian: 



dx 2 



dy 2 



(9) 



The passage from unitless variables (i, y) to variables 
(x, y) with same units modifies the equation only by a 
global factor that docs not affect the global contrast be- 
tween different features on an image plot. 

Similarly to the second derivative, the mean curvature 
function has an equivalent in 2D for a function /(£, y), 
which is given by: 
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When the independent variables carry the same units, 
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we need to use the transformations 

Considering that the spectral function / is defined 
to a factor I , we get: 



[go + (a^) }oy- 
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(ii) 



where a global factor has been removed and Co = (A>£) 2 
is a free positive parameter. 

Let's now compare both 2D methods. In Fig. EJa), 
we plot a Chinese character (hao, which means "good"). 
The character has been broaden by a Gaussian distribu- 
tion and further blurred by a boxcar filter. Although the 
character is recognizable on the raw image, the strokes 
are not sharp. The Laplacian of this image is displayed 
in Fig. EJb). While the Laplacian allows to sharpen the 
strokes a little, the latter remain broad and the whole 
character appears distorted. In contrast, the result ob- 
tained by the 2D curvature method and shown in Fig. 
[31(c) gives a much better representation of the original 
character, with very sharp strokes. Only little distor- 
tion can be observed near stroke intersections and near 
the beginning and the end of each stroke. Analysis of 
real ARPES data with experimental noise leads to simi- 
lar conclusion. In Fig. Eld), we display the ARPES pho- 
toemmission intensity mapping around the Brillouin zone 
center of a Bao.6Ko.4Fe2As2 sample, which has been inte- 
grated over a ± 10 mcV energy range around the Fermi 
level. The high intensity regions represent the Fermi sur- 
face. Although the raw data are sufficient to distinguish 



the presence of two Fermi surface sheets 32j, the Fermi 
surface contours arc difficult to identify precisely. In this 
case, the Laplacian improves the Fermi surface determi- 
nation of the two concentric Fermi surfaces centered at 
the Brillouin zone center. Further improvement is pro- 
vided by the 2D curvature, which makes the Fermi sur- 
face contours narrower. 



B. Inequivalent independent variables 

Unfortunately, spectroscopic data cannot always be 
presented as 2D mappings with x and y axes having 
the same units. This is particularly true when deal- 
ing with the momentum space, like in ARPES, INS and 
RIXS. Commonly, the results may represent the spec- 
tral intensity as a function of energy, and momentum or 
momentum-transfer. In that case, the Laplacian can be 
adapted to variables x and y with different units by using 
the transformations — > and -J^ — » ri-§-, where £ 

dx ^ dx dy ' dy ' ^ 

and rj are positive parameters carrying the same units as 
x and y, respectively. Accounting once more for a global 
positive factor Iq in the absolute value of the experimcn- 
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FIG. 3. Color online, (a) Image representation of the Chinese 
character hao (see the text). (b)[(c)] Corresponding intensity 
plot of the Laplacian [2D curvature]. The original character 
in (a)-(c) is given by the red lines, (d) ARPES Fermi surface 
mapping of Bao.eKo.4Fe2As2. (e)[(f)] Corresponding intensity 
plot of the Laplacian [2D curvature]. 



Where we removed a global factor. The latest equation 
has only one independent parameter, £/rj. A natural 
choice of parameter to capture the main features in an 
image plot is to make the second derivative terms of the 
same order of magnitude, which is done by setting the 
ranges of the data in x and y to similar values. For a 
square grid for example (same number of columns and 
rows), that statement is equivalent to ^/rj = Ax/ Ay, 
where Ax and Ay are the stepsizes along the x and y 
axes, respectively. 

Similarly to the Laplacian, equation (fTU|) can be 
adapted to variables x and y with different units. Us- 
ing the same transformations for x and y, we get: 
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where C x = I 2 ^ 2 and C y = I 2 r] 2 are the only two (pos- 
itive) free parameters for this equation. Using the same 
arguments as for the Laplacian, we can set £/rj = Ax/ Ay 
to assure a good visual representation. In this condition, 
we verify easily that in the limit where Iq — > 0, and thus 
C x —> and C y 0, equation (fT4|) is simplified to 
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The latest equation diverges when: 



(16) 
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(17) 
(18) 



which is equivalent to our definition given in equation 
(|12[) of the Laplacian with variables carrying units. In 
the opposite limit, when Iq oo, we find: 



which corresponds exactly to the position of the extrema 
of /. Therefore, we conclude that the 2D curvature is 
necessarily an improvement compared to the Laplacian 
in tracking the position of extrema. 
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In Figure HI we compare the Laplacian and the 2D cur- 
vature intensity plots for the simulated electronic disper- 
sion given in Figure [TJa). As expected, the 2D curvature 
method gives sharper features. In addition, it tracks the 
original band dispersion with higher accuracy over the 
whole range of energy. It is also instructive to note that 
while the ID curvature method using EDCs and MDCs 
gives results better than the 2D curvature near the band 
bottom and near the Fermi level, respectively, the 2D 
curvature is more reliable over the whole energy range. 




-0.2 0.0 0.2 -0.2 0.0 0.2 

Momentum (nla) Momentum (ji/a) 

FIG. 4. Color online, (a) Laplacian of the simulated ARPES 
intensity plot shown in Figure 1(a). (b) 2D curvature of the 
simulated ARPES intensity plot shown in Figure 1(a). 



broad. Besides, band dispersions are 2D objects (k vs 
E), which are thus more naturally represented by a 2D 
image plot. Indeed, MDC- and EDC-analysis in ARPES 
often lead to slightly different dispersion, even though 
real electronic dispersions, namely E vs k relationships, 
are uniquely defined objects. By using the 2D curva- 
ture method described here, it is possible to remove this 
ambiguity. However, we note that such analysis is accu- 
rate only when we dispose of sufficient data along both 
directions (E and k). 

Although the curvature technique constitutes an ob- 
vious improvement over the second derivative method 
in terms of reliability and sharpness of the spectral fea- 
tures, its main apparent disadvantage is the introduction 
of arbitrary parameters. As shown above, the curvature 
method is at least as reliable as the second derivative 
method in tracking the peak position of dispersive fea- 
tures, whatsoever the parameters used. Similarly, the 
sharpness of the dispersive features is also improved com- 
pared to the second derivative method. In that sense, 
the arbitrariness of the parameters is not a handicap. In 
fact, it gives some latitude to tune the relative contrast 
between different features from a single image and allow 
a better visualization effect. 



V. CONCLUSIONS 



IV. DISCUSSION 

As with the second derivative method, the curvature 
analysis technique described in this paper is a powerful 
method to enhance dispersive features in a spectroscopic 
image. It is very important to keep in mind that this is 
its only purpose and that the information contained in 
the original spectra is indeed richer, despite being some- 
times difficult to access. These visualization methods 
can thus be regarded as effective complementary tools 
in understanding spectroscopic data. For example, while 
the precise shape of MDCs and EDCs from ARPES data 
are often intimately related to intrinsic scattering and 
other electronic interactions, information completely lost 
in the curvature intensity plots, MDCs and EDCs are 
not always good ways to represent dispersion. This is 
especially true for multi-bands systems when bands are 



We have developed a method based on the concept of 
curvature to analyze spectroscopic image plots. As with 
the second derivative method, which is widely used, the 
method presented here is quite efficient for representing 
dispersive features. Using simulated and experimental 
spectral images, we demonstrated that compared to sec- 
ond derivative analysis, the new curvature method im- 
proves significantly the reliability in tracking dispersive 
feature. Moreover, it sharpens spectral features for a 
better visualization of the spectroscopic features. 
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